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Wc present one-dimensional, timc-dcpcndcnt models of the clumps generated by the line- 
deshadowing instability. In order to follow the clumps out to distances of more than 1000 i?*, 
we use an efficient moving-box technique. We show that, within the approximations, the wind 
can remain clumped well into the formation region of the radio continuum. 



1 Introduction 

The line-driven stellar winds of hot stars are subject 
to a strong line-deshadowing instability (e.g. Owocki 
& Rybicki 1984). which causes the wind to become 
highly structured. This structure takes the form of 
strong shocks, strong density contrasts and regions 
of hot, but generally rarefied, gas. 

The structure caused by the line-deshadowing in- 
stability is small-scale and stochastic in nature, as 
opposed to the large-scale, coherent structure asso- 
ciated with discrete absorption components and re- 
lated features in ultraviolet spectral lines of hot stars 
(Prinja 1998) We use the word clumping to refer to 
the small-scale density structure only, with the line- 
deshadowing instability as its most likely cause. 

The degree of clumping at a certain distance r 
from the star is most readily described by the clump- 
ing factor /ci, defined as 
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where <> denotes the time-averaging of the quan- 
tity between brackets. If all of the mass is concen- 
trated in the dense clumps, then the clumping factor 
is the inverse of the volume filling factor, and is equal 
to the overdensity of the clumps with respect to the 
mean wind: 
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The mass-loss rate derived from a density-squared 
dependent observational diagnostic is inversely pro- 
portional to the square root of the clumping factor. 
Most theoretical studies of clumping are limited 
to the wind below 30 stellar radii (i?*). There is, 
however, ample reason to study clumping at much 
larger distances from the star. The radio contin- 
uum used to derive the mass-loss rates of hot stars 
is formed by free-free emission and hence is strongly 
sensitive to clumping. To know the true value of 
the mass-loss rate, we therefore need to know the 
degree of clumping. The same holds for other diag- 
nostics of the mass-loss rate that are proportional 



to the density squared, such as Ha. As is shown 
throughout these proceedings, this is a surprisingly 
difficult thing to do. All mass-loss rate diagnostics 
are affected by uncertainties. One way to reduce 
such uncertainties, is to combine different observa- 
tional mass-loss rate diagnostics, formed in different 
parts of the wind, to obtain the radial stratification 
of clumping. Such a study has been performed by 
Puis et al. (2006). 

Even from such a study, it is impossible to derive 
absolute values of the clumping factor. If one de- 
rives a certain radial stratification of the clumping 
factor assuming the clumping vanishes in the radio 
formation region, then the observations can also be 
explained by this clumping factor multiplied by a 
constant factor, providing the mass-loss rate is low- 
ered accordingly. The derived value of the clumping 
factor (and hence the value of the mass-loss rate) 
thus depends on the assumption one makes about 
the amount of clumping in the radio formation re- 
gion. Therefore it is important to gain insight in how 
clumps evolve as they move out to large distances, 
and to investigate whether clumps can survive as far 
as the radio formation region. 

2 Hydrodynamical models 

2.1 Hydrodynamical models including 
the line-deshadowing instability 

We solve the conservation equations of hydrodynam- 
ics, using the time-dependent hydrodynamics code 
VH-1, developed by J. M. Blondin, and modified by 
S. P. Owocki to include the acceleration due to line 
driving. Our models are one-dimensional. The ra- 
diative acceleration is included in the model using 
the smooth source function method (Owocki 1991). 
The structure is self-excited, in the sense that there 
are no external perturbations at the base of the 
wind. The structure is seeded by internal base per- 
turbations, that arise as radiation is scattered back 
to the base from the structured outer wind. (Initial 
structure arises as the wind solution adapts from the 
smooth initial condition). In the absence of detailed 



M. C. Runacres 



knowledge of the photosphcric perturbations acting 
at the base of a real wind, self-excited structure can 
be seen as a conservative estimate of wind structure. 
Radiative and adiabatic cooling are included in 
the energy equation. Photo-ionisation heating is 
mimicked by imposing a distance-dependent floor 
temperature, below which the temperature is not 
allowed to drop. Details are given in Runacres & 
Owocki (2002). From test calculations performed 
in that paper, we have learned that the amount of 
clumping depends on the value adopted for the floor 
temperature, as well as on the strength of the line- 
driving. Also, it is necessary to maintain a rather 
fine spacing of the radial mesh, in order to ade- 
quately resolve the structure. On the other hand, 
clumping docs not depend on the radiative force be- 
yond 30 i?*. This reduces the outer- wind evolution 
to a pure gasdynamical problem, allowing us to con- 
struct vastly more economical models, which will be 
presented in the next section. 



2.2 Moving-box models 

For a star like C Pup. about half of the radio contin- 
uum is formed beyond 100 i?*. So in order to make 
meaningful predictions about the effect of clump- 
ing on the radio mass loss rate, we need to model 
structure out to very large distances from the star. 
Even without the evaluation of the radiative force, 
evolving the entire stellar wind (between 1 and say 
1000 i?*) at the required high spatial resolution is 
still very expensive. A solution is suggested by re- 
alising that the structure generated by the instabil- 
ity, apart from being stochastic, is also quasi-regular 
in the sense that similar features are repeated over 
time. Therefore it is not necessary to keep track of 
the whole stellar wind during the duration of the 
simulation. It is enough to select a limited but rep- 
resentative portion of the structure, and follow this 
"box" as it moves out at the terminal speed. Follow- 
ing a portion of the wind entails transforming the 
conservation equations to a moving reference frame. 
This is not possible directly, as the spherical equa- 
tions of hydrodynamics are not invariant under a 
Galilean transformation. This problem can be cir- 
cumvented by rewriting the equations in a pseudo- 
planar form. In this form, the equations resemble 
the planar equations of hydrodynamics, while still 
describing a spherical geometry. We impose peri- 
odic boundary conditions on the box, i.e. structures 
that flow out of the box on one side, are made to 
enter it on the other side. Details can be found in 
Runacres & Owocki (2005). 

In the following section, we use a periodic box 
model, starting from a hydrodynamical model in- 
cluding the line-dcshadowing instability, to predict 
the radial stratification of wind clumping. The 
adopted model parameters are the same as in 
Runacres & Owocki (2005). 



3 Results 



Fig. 1 shows the density contrast (density divided 
by mean density) within the box as a function of 
radius and time, as the box moves out from ~ 100 
to r^ 1300 i?*. The backward running streaks are 
shells that are somewhat slower than the terminal 
speed, the forward running streaks are faster than 
the terminal speed. The streaks broaden as they 
evolve, reflecting the fact that shells expand (at a few 
times the sound speed) as they move out. Within the 
assumptions of the model, the dumpiness is main- 
tained by collisions between shells. As shells collide, 
they form denser shells, counteracting their pressure 
expansion. 
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Figure 1: Image of the density contrast as a function 
of space and time. The intensity scale has 
been truncated to highlight the kinematics 
of the shells. The adopted model parame- 
ters are the same as in Runacres & Owocki 
(2005). 



The clumping factor for a moving-box model ex- 
tending out to ~ 1300 i?* is shown in Fig. 2. Be- 
low 100 i?* the model is a line-driven instability 
model, above 100 i?* it is a moving-box model. It is 
clear that these models predict that the winds stays 
clumped well into the radio formation region, with 
clumping factors beyond 200 i?* ranging from 2.5 to 
6. Inferred mass- loss rates would therefore be over- 
estimated by a factor of two. 
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Figure 2: The clumping factor as a function of radius, predicted by our theoretical models. Below 100 i?* 
the model is a line-driven instability model, above 100 i?» it is a moving-box model. 



4 Discussion and conclusions 

Our models predict an increase of the clumping fac- 
tor from the base of the wind to ^ 50 R^, (Fig. 2), 
after which the clumping factor decreases, maintain- 
ing a level of residual clumping beyond 200 i?* . This 
does not quite match the clumping factors derived 
from observations by Puis et al. (2006), which start 
to tail off closer to the star (~ 10 i?*). As has been 
mentioned before, the observations do not tell us 
whether or not there is residual clumping at very 
large distances from the star. 

There are of course a number of limitations to 
our model. As has been mentioned above, we have 
used self-excited structure without external pertur- 
bations. Also, we have not attempted to model dif- 
ferent spectral types. In particular, the important 
difference between dense and less dense winds found 
by Puis et al. (2006) has not been investigated at 
all in our models. 

A key limitation of the present model is of course 
its restriction to just one dimension. The focus here 
is entirely on the extensive radial structure, and the 
instabilities that are likely to break up the azimuthal 
coherence of the structure are not accounted for. It 
remains to be seen to what extent this changes the 



global evolution of instability-generated wind struc- 
ture. We plan to extend the present models to 2-D 
in the near future. 
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